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In  order  to  investigate  the  effect  of  capillary  pressure  on  the  transport  of  liquid  water  in  the  cathode 
gas  diffusion  layer  (GDL)  of  a  polymer  electrolyte  fuel  cell,  a  one-dimensional  steady-state  mathematical 
model  was  developed,  including  the  effect  of  temperature  on  the  capillary  pressure.  Numerical  results 
indicate  that  the  liquid  water  saturation  significantly  increases  with  increases  in  the  operating  temper¬ 
ature  of  the  fuel  cell.  An  elevated  operating  temperature  has  an  undesirable  influence  on  the  removal  of 
liquid  water  inside  the  GDL.  A  reported  peculiar  phenomenon  in  which  the  flooding  of  the  fuel  cell  under 
a  high  operating  temperature  and  an  over-saturated  environment  is  more  serious  in  a  GDL  combined  with 
a  micro-porous  layer  (MPL)  than  in  a  GDL  without  an  MPL  [Lim  and  Wang,  Electrochim.  Acta  49  (2004), 
4149-4156]  is  explained  based  on  the  present  analysis. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

In  polymer  electrolyte  fuel  cells  (PEFC),  excessive  liquid  water 
remaining  inside  the  gas  diffusion  layer  (GDL),  catalyst  layer  (CL) 
or  even  in  the  gas  channel  tends  to  inhibit  the  performance  of  the 
fuel  cell,  especially  at  a  high  current  density.  This  phenomenon, 
commonly  called  flooding  [1],  is  confirmed  to  be  more  serious  in 
the  cathode  GDL  than  in  the  GDL  on  the  anode  side  [2].  Though  the 
GDL  is  commonly  wet-proof  treated  by  adding  polytetrafluroethy- 
lene  (PTFE)  to  increase  its  hydrophobicity  or  combined  with  a  very 
highly  hydrophobic  micro-porous  layer  (MPL),  the  flooding  still 
appears  or  is  even  more  serious  under  some  performance  condi¬ 
tions  than  it  with  a  lower  PTFE  content  or  in  a  GDL  without  an  MPL 
[3,4].  Thus,  further  understanding  of  the  mechanism  of  liquid  water 
transport  inside  the  GDL  is  very  critical. 

Within  the  past  decade,  several  experimental,  numerical,  and 
analytical  investigations  of  water  removal  have  been  performed 
[5-9].  There  are  two  kinds  of  opinions  regarding  the  mechanism 
of  liquid  water  transport  inside  the  hydrophobic  GDL.  Nam  and 
Kaviany  suggested  based  on  their  theoretical  analysis  that  the  dis¬ 
tribution  of  water  showed  a  branching- type  geometry  [5].  They 
considered  that  water  vapor  condenses  at  the  surface  into  micro¬ 
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droplets,  which  intermittently  agglomerate  to  macro-droplets  of 
liquid  water.  In  their  model,  the  large  main  streams  extend  and 
branch  into  many  small  streams  from  the  CL  to  the  channel  side.  A 
similar  model  was  also  described  by  Pasaogullari  and  Wang  [7].  This 
kind  of  model  was  termed  an  “upside-down  tree”  or  “inverted  tree” 
capillary  network  by  Litster  and  Djilali  [9,10],  in  order  to  distinguish 
it  from  another  transport  model  that  they  proposed  themselves 
based  on  their  ex  situ  experimental  observation  using  a  fluores¬ 
cence  microscope  [9].  In  contrast  to  the  “upside-down  tree”  model, 
they  proposed  a  mechanism  in  which  the  accumulation  of  liquid 
water  generated  by  the  oxygen  reduction  reaction  (ORR)  builds  up 
liquid  pressure  and  preferentially  breaks  through  the  pores  with  the 
greatest  diameter.  A  dominated  pathway  is  formed  after  the  liquid 
water  penetrates  through  the  GDL  along  the  thickness  direction. 
Meanwhile,  numerous  “dead  ends”  exist  inside  the  GDL,  where  an 
adjacent  breakthrough  channel  forms.  In  fact,  a  similar  model  was 
proposed  by  Yang  et  al.  [11  ]  earlier  than  Litster  et  al.  [9].  Yang  et  al. 
believed  that  some  preferential  openings  are  connected  to  a  net¬ 
work  of  large  hydrophilic  pores  inside  the  GDL  to  form  a  pathway 
for  liquid  water  transport  [11  ]. 

Though  these  two  models  differ  in  their  details,  they  share  the 
concept  that  the  capillary  action  resulting  from  the  capillary  pres¬ 
sure  distribution  plays  a  dominant  role  in  the  transport  of  liquid 
water  inside  the  GDL.  Djilali  [10]  also  proposed  that  liquid  water 
transport  is  dominated  by  means  of  capillary  diffusion  because 
the  capillary  number  (the  ratio  of  viscous  forces  to  surface  tension 
forces)  is  generally  small. 


0378-7753 /$  -  see  front  matter  ©  2008  Elsevier  B.V.  All  rights  reserved. 
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Nomenclature 

F 

Faraday  constant  (C  mol-1 ) 

h 

heat  transfer  coefficient  (W(m2 1<)-1 ) 

I 

current  density  (A cm-2) 

j 

Leverett  function 

I< 

permeability  (m2) 

I< rl 

relative  permeability 

Mh20 

molecular  weight  of  water  (kg  mol-1 ) 

He- 

stoichiometry  coefficient  of  electron  in 

cathode 

reaction 

Nu 

Nusselt  number 

Pc 

capillary  pressure  (Pa) 

Q 

heat  flux  (Wm-2) 

s 

saturation  of  liquid  water 

ASC 

entropy  change  of  oxygen  reduction 

reaction 

(J(molK)-1) 

T 

temperature 

X 

x  coordinate  (m) 

Greek  symbols 

P 

constant  in  capillary  pressure  versus  temperature 

equation  (I<-1) 

k 

thermal  conductivity  (W  nrr 1  K-1 ) 

8 

thickness  (m) 

£ 

porosity  of  porous  medium 

over-potential  (V) 

0 

contact  angle 

V 

kinematic  viscosity  (m2  s-1 ) 

o 

surface  tension  (N  m-1 ) 

Subscripts 

c 

cathode 

f 

fluid 

g 

gas  phase 

1 

liquid  phase 

L 

land  of  bipolar  plate,  or  left  hand  side 

MPL 

micro-porous  layer 

nw 

non-wetting  phase 

r 

reference  value 

R 

right  hand  side 

Based  on  this  common  idea,  some  investigations  of  the  liq¬ 
uid  water  distribution  in  GDL  and  the  effect  of  capillary  pressure 
were  carried  out  numerically  or  analytically.  Nam  and  Kaviany  [5] 
assumed  that  water  vapor  (or  steam)  is  produced  by  a  reaction  on 
the  cathode  CL  surface.  The  phase  change  of  the  water  accompa¬ 
nies  its  transport.  Driven  by  capillary  pressure,  condensed  liquid 
water  flows  from  a  location  with  higher  liquid  water  saturation 
(volume  fraction  of  the  open  pore  spaces  occupied  by  liquid  water) 
to  a  point  with  lower  saturation.  They  analyzed  the  liquid  water 
saturation  distribution  and  the  effects  of  the  porosity  and  per¬ 
meability  of  the  GDL.  They  also  compared  the  capillary  pressure 
and  liquid  water  saturation  in  a  GDL  with  an  MPL  (micro-porous 
layer)  with  those  in  a  GDL  without  an  MPL.  In  contrast,  Pasaogullari 
and  coworkers  [2,7]  considered  that  liquid  phase  water  was  pro¬ 
duced  on  the  cathode  CL  surface.  The  phase  change  was  neglected 
in  their  one-dimensional  model.  They  obtained  the  liquid  water 
saturation  distribution  inside  the  GDL  and  MPL  by  an  analytical 
solution  and  discussed  the  effect  of  the  contact  angle  of  the  water. 
Similar  work  was  performed  by  Zhan  et  al.  [12,13].  These  analy¬ 
ses  were  based  on  a  simplified  one-dimensional  model;  however, 
they  can  predict  the  liquid  water  saturation  distribution  inside  the 


GDL  along  the  thickness  direction,  which  is  very  difficult  up  to  now 
either  by  optical  diagnostics  [9,11,14,15]  or  by  neutron  radiography 
measurement  [16-18].  By  two-dimensional  numerical  simulation, 
Natarajan  and  Nguyen  [19]  calculated  the  liquid  water  saturation 
distribution  using  a  capillary  pressure-saturation  relationship  pro¬ 
posed  by  them.  They  found  that  the  saturation  at  the  catalyst  layer 
is  close  to  100%.  But  they  neglected  the  effect  of  temperature  on 
capillary  pressure.  Berning  and  Djilali  [20]  found  that  the  liquid 
water  saturation  decreases  from  the  catalyst  layer  towards  the  gas 
channel,  which  is  similar  to  that  reported  by  You  and  Liu  [21  ]. 

In  a  porous  medium,  the  capillary  pressure  is  related  to  the 
surface  tension  and  the  contact  angle  (or  wetting  coefficient)  of 
the  fluid  on  the  solid  surface.  It  is  important  to  emphasize  that 
both  surface  tension  and  contact  angle  are  decreasing  functions  of 
temperature  [22-24].  Consequently,  the  capillary  pressure  in  the 
porous  medium  tends  to  decrease  with  increases  in  temperature 
[22-24].  For  example,  the  relative  decrease  in  capillary  pressure  is 
typically  about  -1%K-1  in  porous  media  samples  fashioned  from 
soils,  rocks,  etc.  [23].  In  most  of  the  reported  analytical  solutions 
or  two-phase  numerical  models  of  liquid  water  transport  in  the 
PEFC,  only  the  effect  of  temperature  on  the  surface  tension  of 
the  air-water  system  is  usually  considered  [25,26],  whereas  the 
influence  of  temperature  on  the  contact  angle  is  usually  ignored. 
Recently,  Kumbur  et  al.  [27]  experimentally  measured  the  capillary 
pressure-saturation  relationship  for  several  commercial  GDL  sam¬ 
ples  at  different  temperatures  and  found  that  the  capillary  pressure 
decreases  with  increases  in  temperature.  In  an  empirical  corre¬ 
lation  of  the  capillary  pressure-saturation  relationship  that  they 
proposed  the  effect  of  temperature  was  introduced  explicitly.  One 
distinctive  feature  of  their  empirical  correlation  is  that  no  con¬ 
tact  angle  parameter  is  required  as  an  input  except  for  the  PTFE 
weight  percentage  of  GDL.  They  claimed  that  their  empirical  corre¬ 
lation  successfully  predicted  the  capillary  pressure  over  the  entire 
range  of  saturation.  Flowever,  the  figures  in  their  paper  show  that 
their  correlation  is  not  highly  consistent  with  their  experimental 
data  for  most  of  the  GDL  samples  used  for  their  experiments  when 
the  saturation  is  less  than  0.1.  In  particular,  when  the  saturation  is 
very  small,  such  as  less  than  0.01,  their  empirical  correlation  pre¬ 
dicts  negative  capillary  pressure  values  for  a  hydrophobic  GDL  with 
5  wt.%  PTFE,  which  is  physically  incorrect.  In  fact,  for  a  GDL  working 
in  an  actual  PEFC,  the  water  saturation  is  mostly  less  than  0.1,  and 
serious  flooding  occurs  if  the  saturation  is  higher  than  0.1  [7].  In 
contrast,  the  standard  Leverett  function  can  avoid  such  error  even 
for  very  low  levels  of  saturation,  although  it  is  originally  derived 
from  soil  science. 

In  this  work,  we  propose  a  one-dimensional  steady-state 
mathematical  model  based  on  a  capillary  pressure-saturation  rela¬ 
tionship,  including  the  influence  of  temperature,  using  the  standard 
Leverett  function.  The  influence  of  temperature  on  the  liquid  water 
saturation  is  investigated  numerically  inside  the  cathode  GDL.  The 
phenomenon  related  to  liquid  water  transport  reported  in  the  lit¬ 
erature  is  discussed  based  on  the  present  analyses. 


2.  Physical  and  mathematical  models 

As  mentioned  in  the  previous  section,  flooding  usually  occurs 
more  seriously  in  the  GDL  at  the  cathode  side  than  that  at  the  anode 
side.  Thus,  the  present  study  only  considers  the  water  transport  in 
the  GDL  at  the  cathode  side.  A  schematic  physical  model  is  shown  in 
Fig.  1 .  In  the  cathode  catalyst  layer,  the  water  and  heat  are  generated 
due  to  the  oxygen  reduction  reaction: 

4H+  +  02(g)  +  4e_  2H20(1). 


(i) 
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Fig.  1.  Physical  model. 


.Channel 


Here,  the  generated  water  is  assumed  to  be  in  the  liquid  phase. 
The  phase  change  of  the  water  is  neglected  for  the  purpose  of  sim¬ 
plification;  thus,  the  water  substance  is  transported  in  the  GDL 
in  the  liquid  single  phase.  Heat  convection  due  to  gas  and  liquid 
water  flow  inside  the  GDL  is  also  ignored.  The  CL  is  assumed  to  be 
an  infinitely  thin  layer,  and  is  also  assumed  to  act  as  an  adiabatic 
boundary  because  of  its  small  thermal  conductivity  in  comparison 
with  that  of  the  GDL  [28].  Therefore,  all  of  the  heat  generated  at 
the  catalyst  layer  is  transferred  to  the  gas  flowing  in  the  gas  chan¬ 
nel  and  bipolar  plate  by  heat  conduction.  Under  a  steady  state,  a 
one-dimensional  energy  conservation  equation  and  boundary  con¬ 
ditions  are  expressed  as  follows: 


At  x  =  0  : 

At  x  =  8  : 


(2a) 

f-AScTo 

q={  ne-F  +  *)L 

(2b) 

h(T$  -  Tf)  XL(TS-T[) 
q  2  2<$l 

(2c) 

In  Eq.  (2b),  A Sc  is  the  entropy  change  in  the  reaction 
of  Eq.  (1),  which  represents  the  reversible  heat  generation. 
ASC  =  -326.36 J  (mol  K)-1  indicates  that  the  reaction  is  exothermic 
[29].  ?7c  is  the  over-potential  of  the  cathode  CL,  which  represents  the 
heat  generation  due  to  irreversible  operation  (i.e.,  effective  losses) 
[28,30].  Where  T0  is  the  temperature  at  the  interface  between  the  CL 
and  GDL  (i.e.,  at  x  =  0),  which  equals  the  temperature  of  the  CL.  ne- 
is  the  stoichiometry  coefficient  of  the  electron  and  ne-  =  4  from  Eq. 
(1).  /  is  the  average  current  density.  On  the  GDL  surface  facing  the 
gas  (air)  channel  (i.e.,  at  x  =  8,  the  heat  convection  is  represented  by 
heat-transfer  coefficient  h,  which  is  given  from  the  Nusselt  number 
(Nu)).  The  value  of  h  is  calculated  to  be  106.1  W  (m2  K)-1  for  the  air 
flow  because  Nu  is  constant  and  equals  3.61  for  laminar  flow  in  the 
square  duct  under  constant  heat  flux  [31  ].  Ts  and  Tf  are  the  tempera¬ 
tures  at  the  GDL-channel  interface  (x  =  8)  and  the  bulk  temperature 
of  the  gas  in  the  gas  channel,  respectively.  The  temperature  at  the 
right  side  wall  of  the  bipolar  plate  is  also  assumed  to  be  Tf.  The  ther¬ 
mal  conductivities  of  the  GDL,  bipolar  plate  and  air  in  the  channel 
are  listed  in  Table  1 .  For  a  porous  GDL,  an  effective  thermal  conduc¬ 
tivity  of  the  pores  and  solid  regions  in  combination  is  applied  [37]. 
Here,  a  low  thermal  conductivity  of  the  GDL  is  utilized,  considering 
the  fact  that  the  thermal  conductivity  of  the  GDL  decreases  with 
increases  in  temperature  and  PTFE  content  [38]. 


Table  1 

Parameters  of  calculation 


Geometry  size  (mm) 

Thickness  of  GDL  with/without  MPL  (5) 

0.27/0.3 

Channel  width/height 

1/1 

Bipolar  plate  thickness  <5L/height 

1.5/2 

Thickness  of  MPL 

0.03  [2] 

Thermal  conductivity  (W  m- 1  K-1 ) 

GDL  (A.)  (effective) 

0.75  [32] 

Air  (Aair) 

0.0294  [33] 

Bipolar  plate  (AL) 

20  [26] 

Faraday  constant  (F) 

96487  C  moU1 

Molecular  weight  of  water  (MH2 0 ) 

18.0  x  10-3  kg  mol-1 

Surface  tension  of  liquid  water  at  25  °C  (cr) 

0.0725  Nm-1  [34] 

Contact  angle  of  water  at  25  °C  unless 
specifically  mentioned  ( 6 ) 

112° 

Viscosity  of  liquid  water  (v) 

3.5  x  10-7  m2  s-1  [2] 

Porosity  of  GDL  (e) 

0.78  [34] 

Porosity  of  MPL  (£MPL) 

0.5  [2] 

Single  phase  permeability  of  GDL  (K) 

8.99  x  10-12  m2  [35] 

Single  phase  permeability  of  MPL  ( I<mpl  ) 

2.47  x  10“16  m2  [2] 

Gas  temperature  in  channel  (Tf) 

80  °C 

Cathode  over-potential  (rjc) 

0.4  V  [36] 

Current  density  (/) 

1.0  A  cm-2 

The  capillary  pressure  is  defined  as  the  difference  between  liq¬ 
uid  and  gas-phase  pressures: 


Pc  =  Pi-Pg. 


(3) 


If  the  gas  phase  pressure  inside  the  GDL  is  assumed  to  be  con¬ 
stant  at  Pg,  which  is  equal  to  the  environmental  pressure,  Pc  may  be 
considered  a  relative  pressure  head  of  liquid  water.  In  hydropho¬ 
bic  porous  media,  Pc  >  0,  i.e.,  the  liquid-phase  pressure,  is  higher 
than  the  gas-phase  pressure,  whereas  in  hydrophilic  media,  the 
liquid-phase  pressure  is  lower  than  the  gas-phase  pressure  (Pc  <  0). 
Furthermore,  in  hydrophobic  porous  media,  the  magnitude  of  the 
capillary  pressure  increases  with  the  increase  of  the  liquid  water 
saturation.  Therefore,  a  liquid  pressure  drop  is  formed  from  high  to 
low  liquid  saturation  regions,  which  is  an  important  driven  force 
for  liquid  water  transport.  In  many  reports  [2,7,25,26],  Eq.  (4)  is 
adopted  for  the  capillary  pressure  as  a  function  of  non-wetting 
phase  saturation  (snw)- 


Pc(Tr)  =  ~  gC°S^2J(Snw).  (4) 

(K/e)1/2 

The  minus  symbol  is  imposed  to  match  with  our  definition  of 
capillary  pressure,  i.e.,  Eq.  (3).  In  Eq.  (4)  J(snw)  is  the  Leverett  func¬ 
tion,  which  was  originally  regressed  by  Udell  [39]  based  on  the 
experimentally  obtained  capillary  pressure  data  sets  in  unconsol¬ 
idated  sands  reported  by  Leverett  [40].  The  Leverett  function  was 
extended  to  hydrophobic  porous  media  by  Nam  and  Kaviany  [5] 
and  Pasaogullari  and  Wang  [7]  for  the  fuel  cell  transport  model.  The 
snw  is  defined  as  the  volume  fraction  of  the  void  occupied  by  the 
non-wetting  phase.  In  hydrophobic  media,  the  non-wetting  phase 
is  liquid  water,  i.e.,  snw  =  s,  whereas  in  hydrophilic  media,  the  non¬ 
wetting  phase  is  air,  and  snw  =  1  -  s.  Here,  s  indicates  the  saturation 
of  liquid  water,  which  is  defined  as 


^liquid 

^pore 

Then,  the  J(s)  function  is  expressed  as  [7]: 


(5) 


J(s)  =  < 


1.417(1  —  s)  —  2.12(1  -s)2 
+1.263(1  -s)3, 

1.417s- 2.12s2  +  1.263s3, 


for  0  <  0  <  90° 
for  90°  <6  <  180° 


(6) 
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where  0  is  the  contact  angle  of  liquid  on  a  flat  surface.  For 
hydrophilic  surfaces,  0  ranges  from  0<$<90°,  while  90°  <0<  180° 
for  hydrophobic  surfaces. 

In  the  electrolyte  membrane,  the  water  flux  driven  by  electro- 
osmotic  drag  is  assumed  to  be  balanced  by  backing  diffusion  at 
steady  state,  i.e.,  there  is  no  water  leakage  from  the  electrolyte 
membrane.  All  of  the  water  generated  in  CL  must  be  transported 
through  GDL  by  capillary  pressure  toward  the  gas  channel.  For  sim¬ 
plification,  assuming  that  the  phase  change  of  water  inside  the 
GDL  and  the  effect  of  gravitation  are  negligible,  Darcy’s  law  can 
be  expressed  as  [7]: 

^Mh2o  =  -^VPc(T),  (7) 

where  F  is  the  Faraday  constant,  MH2o  is  the  molecular  weight  of 
water,  v  is  the  viscosity  of  liquid  water,  I<  is  the  absolute  perme¬ 
ability  of  GDL  and  kT\  is  the  relative  permeability  for  liquid  water. 
The  relative  permeability  can  be  determined  by  means  of  experi¬ 
ment  [34]  or  pore  network  modeling  [41,42].  In  this  work,  the  GDL 
is  assumed  to  be  isotropic  and  homogeneous  porous  medium,  and 
the  relative  permeability  of  individual  phase  is  assumed  to  be  pro¬ 
portional  to  the  cube  of  individual  phase  saturation  [7].  For  liquid 
water,  it  is  known  as 

krl=S3.  (8) 

A  boundary  condition  for  Eq.  (7)  is  given  at  the  GDL  surface  facing 
the  gas  channel  (at  x  =  8).  In  general,  a  constant  value  is  assigned  for 
the  liquid  water  saturation  atx  =  <$  [5,7,12,13].  Here,s  =  0  is  assumed 
at  x  =  8: 

at  x  =  8  :  s  =  0.  (9) 

The  effect  of  temperature  on  capillary  pressure  is  expressed  by 
a  formula  proposed  by  Grant  [22,23],  as  follows: 

PcCn  =  PcCrr)!±L  (10) 

where  Pc(Tr)  is  capillary  pressure  at  a  reference  temperature  (TV). 
Here,  Tr  =  298  K  is  used  for  the  following  analysis,  p  is  a  constant 
related  to  porous  medium  [22].  Due  to  lack  of  experimental  data 
for  GDL,  the  constant  p  =  -380.4  K  was  adopted,  which  was  deter¬ 
mined  from  experimental  data  of  air-water  drainage  in  a  silt  loam 
[22].  Fig.  2  shows  the  Pc  versus  temperature  profile.  It  is  worthwhile 
to  emphasize  that  the  influence  of  temperature  on  both  surface 
tension  and  contact  angle  have  been  considered  in  Eq.  (10)  [22]. 
Therefore,  it  is  not  necessary  to  consider  the  effect  of  temperature 
on  the  surface  tension  or  the  contact  angle  individually,  which  is  an 
important  advantage  because  the  effect  of  temperature  on  the  con¬ 
tact  angle  of  GDL  is  not  well  known.  Recently,  similar  method  has 
been  adopted  by  Kumbur  et  al.  [27].  Thus,  Eq.  (10)  is  employed  pri- 
orly  in  the  following  calculations,  i.e.,  the  constant  surface  tension 


and  the  constant  contact  angle  values  at  the  reference  temperature 
(298  K)  are  required  in  the  Eq.  (4)  unless  specified  otherwise.  Eq. 
(10)  is  introduced  into  Eq.  (7),  and  then  liquid  water  conservation 
is  expressed  as 


/  n  n  K\K 

_Mh2o  =  -  — v 


P  +  T 

P  +  T, 


Pc(Tr) 


(11) 


In  the  analysis  of  GDL  combined  with  MPL,  the  liquid  water  sat¬ 
uration  at  the  interface  between  MPL  and  GDL  is  determined  from 
phenomena  of  capillary  equilibrium  expressed  by  the  following 
equation: 


cos(P) 

(K/e)1/2 


J(s)  = 


cos(Pmpl) 

(^mpl/^mpl)1^2 


-J(smpl)- 


(12) 


3.  Results  and  discussion 


3.1.  Temperature  field  in  GDL 
The  analytical  solution  of  (2)  is 

T  =  -lx  +  T0.  (13a) 

Tf  +  [5/A.  +  1/(0. 5h  +  0.5A.L/<5L)]/7cl 

0  1  +  p  +  l/(0.5/i  +  0.51L/4)]ASc//(ne-F)'  v  1 

Under  the  operation  condition  given  by  Tf  =  353K,  rjc  =  0.4  V, 
1=  1.0  A  cm-2  and  the  other  parameters  listed  in  Table  1,  we  obtain 
T0  =  357.0  K  and  Ts  =  354.2  K,  and  the  temperature  field  in  GDL  is 
shown  in  Fig.  3.  The  temperature  difference  between  the  two  sides 
of  the  GDL  is  2.8  K,  which  falls  in  the  same  range  of  temperature 
drop  as  that  described  in  the  literature  [38]. 


3.2.  Saturation  and  capillary  pressure  inside  GDL 

Substituting  Eqs.  (4)  and  (6)  into  Eq.  (11),  the  liquid  water  satu¬ 
ration  in  a  one-dimensional  model  is  expressed  as  follows: 

s3(— 3.789s2  +  3.338s  -  0.966)  ^  =/,  for  6  <  90°  (14a) 

s3(3.789s2  -  4.24s  +  1 .417)  ^  =/,  for  9  >  90°,  (14b) 

where 


A(P  +  Tr)  J(s)s 3  q 
1  P  +  T  +( P  +  T)X ’ 

and 

^  _  JMh2o  v 

2f  {Ks^^cx  cos(0) 


(14c) 


(14d) 
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Fig.  4.  (a)  Liquid  water  saturation  and  (b)  capillary  pressure  profiles  across  the  GDL 
thickness  at  298 1<  and  353 1<.  Solid  lines:  the  temperature  gradient  is  considered 
in  Eq.  (14c).  Dash-dot  lines:  the  temperature  gradient  is  neglected,  i.e.,  q  =  0  in  Eq. 
(14c). 


0  100  200  300 


x  [pm] 

Fig.  5.  (a)  Liquid  water  saturation  and  (b)  capillary  pressure  profiles  across  the 
thickness  of  the  GDL  with  an  MPL  at  298  K  and  353  K. 


The  local  saturation,  s(x),  can  be  obtained  by  numerically  solv¬ 
ing  (14).  The  capillary  pressure  distribution  can  then  be  calculated 
from  Eq.  (10).  For  example,  when  i/c  =  0.4V,  /=1.0Acm-2  and  the 
contact  angle  0  =  112°,  the  liquid  water  saturations  under  different 
temperature  conditions  are  as  shown  in  Fig.  4a.  Since  the  capillary 
pressure  decreases  with  an  increase  in  temperature,  s  at  Tf  =  353  K 
(air  temperature  in  the  gas  channel)  is  significantly  higher  than 
that  at  the  room  temperature  of  air  (Tf  =  298  K).  In  particular,  s(0), 
i.e.,  the  water  saturation  at  the  interface  between  CL  and  GDL  (at 
x  =  0),  for  Tf  =  353 1<  is  about  41.9%  higher  than  that  for  Tf  =  298  K.  The 
capillary  pressure  profile  changes  as  shown  in  Fig.  4b.  At  the  inter¬ 
face  between  CL  and  GDL  (at  x  =  0),  Pc(0)  for  Tf  =  353  K  is  as  low  as 
41.0%  of  that  for  Tf  =  298  K.  This  result  clearly  indicates  that  eleva¬ 
tion  of  the  operating  temperature  significantly  reduces  the  liquid 
water  removal  ability  of  porous  GDL,  increases  the  water  satura¬ 
tion  in  the  entire  volume  of  GDL,  suppresses  the  gas  permeability 
and  eventually  causes  flooding.  All  of  these  influences  decrease  fuel 
cell  performance.  This  effect  of  operating  temperature  has  not  been 
well  recognized  in  the  reported  literature  [2,5,7,12,13],  although 
Weber  and  Newman  found  that  the  performance  of  fuel  cell  may 
become  worse  at  high  operating  temperature  [28]. 

By  solving  (14)  without  the  temperature  gradient,  i.e.,  neglect¬ 
ing  the  second  term  in  Eq.  (14c),  we  can  evaluate  the  effect  of  the 
non-uniformity  of  temperature  on  the  distributions  of  water  sat¬ 
uration  and  capillary  pressure  as  shown  by  the  dash-dot  lines  in 
Fig.  4.  The  results  indicate  that  the  temperature  gradient  in  the  GDL 
has  a  small  influence  on  the  saturation  and  capillary  pressure  dis¬ 
tributions  because  of  the  small  temperature  difference  despite  the 
large  temperature  gradient  in  the  GDL  (9.3  x  103  Km-1 ).  We  spec¬ 
ulate  that  the  similar  conclusion  is  also  valid  for  anode  GDL,  MPL, 


membrane  as  well  as  catalyst  layer  because  of  their  thin  thickness 
too.  Therefore,  we  neglect  the  effect  of  temperature  gradient  in  the 
following  calculation. 

It  is  well  known  that  a  micro-porous  layer  with  high  hydropho- 
bicity,  low  permeability  and  small  pore  diameter  can  usually 
significantly  improve  the  limitation  of  mass  loss  and  effectively 
avoid  flooding  [43,44].  As  an  example,  an  MPL  with  a  thickness 
<Wl  =  30  |jim  is  introduced  here,  and  the  operating  temperature  of 
the  fuel  cell  is  assumed  to  be  uniform  at  298  K.  The  contact  angle 
0  =  112°  for  the  GDL  and  0  =  120°  for  the  MPL,  and  the  surface  ten¬ 
sion  a  =  0.0725  N  nrr 1  at  298  K.  The  other  parameters  are  tabulated 
in  Table  1 .  At  the  interface  between  the  MPL  and  GDL,  the  satu¬ 
ration  must  satisfy  Eq.  (12).  As  shown  in  Fig.  5a,  the  liquid  water 
saturation  inside  the  main  part  of  the  MPL  is  higher  than  that  in 
the  GDL.  A  very  large  capillary  pressure  gradient  is  generated  near 
the  interface  between  the  MPL  and  GDL,  as  shown  in  Fig.  5b,  which 
is  effective  for  draining  liquid  water  from  the  CL  through  the  MPL 
and  toward  the  gas  channel.  The  high  capillary  pressure  is  one  of 
the  significant  advantages  in  comparing  with  other  factors,  such  as 
low  ohmic  losses  [44]. 

As  a  typical  example  referred  to  in  some  reports  [2,43,44],  Qi 
and  Kaufman  [43]  found  that  MPL  significantly  improves  the  per¬ 
formance  of  fuel  cells  at  a  low  cell  operating  temperature  (308  K), 
increasing  the  limiting  current  up  to  0.6  A  cm-2  at  0.4  V  when  air  is 
utilized  as  a  cathode  reactant.  In  contrast,  Lim  and  Wang  [3]  found 
that,  in  a  peculiar  phenomenon,  the  MPL  reduces  the  performance 
of  the  fuel  cell,  i.e.,  the  limiting  current  density  of  the  fuel  cell 
is  lower  (about  0.33  A  cm-2  at  0.4  V)  than  that  without  the  MPL, 
although  the  fuel  cell  was  operated  at  similar  reactant  stoichiom¬ 
etry  and  air-fed  as  the  cathode  reactant,  using  similar  wet-proof 
carbon  paper  as  the  GDL  except  at  the  higher  operating  tempera¬ 
ture  (353 1<)  and  in  an  over- saturated  environment  (hydrogen  and 
air  humidified  at  368  K  and  363  K,  respectively).  Such  a  low  limita- 
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tion  current  indicates  that  flooding  occurred.  They  explained  briefly 
that  the  MPL  imposed  an  additional  diffusion  resistance  to  oxygen 
transport  into  the  catalyst  layer,  and  the  mass  transport  rate  lim¬ 
ited  the  current  [3  ].  However,  they  did  not  explain  the  reason  for  the 
additional  diffusion  resistance.  The  higher  current  density  under  a 
low  operating  temperature  indicates  that  the  reaction  proceeded 
faster  and  that  more  water  was  generated  than  in  the  latter  case 
operated  at  a  higher  temperature.  However,  flooding  occurred  at 
the  latter  case  regardless  of  its  lower  water  generation  rate.  This 
may  be  attributed  to  the  decrease  of  the  driven  force  of  liquid  water 
transportation  in  the  MPL  and  GDL,  i.e.,  the  capillary  pressure  gradi¬ 
ent  at  the  higher  temperature.  Under  a  high  operating  temperature, 
even  under  a  low  reaction  rate,  the  generated  water  could  not  be 
drained,  but  rather  accumulated  inside  the  GDL  and  MPL,  which 
increased  the  diffusion  resistance  of  oxygen  transport  and  further 
decreased  the  reaction  rate. 

In  fact,  in  Lim  and  Wang’s  work  [3],  it  was  found  that  the  con¬ 
tact  angle  of  GDL  rapidly  decreased  from  98°  to  108°  (hydrophobic) 
to  80°  (hydrophilic)  as  the  temperature  increased  from  298 1<  to 
353  K.  At  temperatures  above  313  K,  even  a  40-wt.%  FEP-treated 
hydrophobic  GDL  became  hydrophilic  [3].  Although  their  measure¬ 
ment  of  the  contact  angle  was  focused  only  on  a  normal  GDL,  here, 
we  assume  that  a  similar  phenomenon  occurs  in  the  MPL.  Recently, 
Kumbur  et  al.  [27]  experimentally  found  that  the  hydrophilic  pore 
volume  of  hydrophobic  treated  GDL  increased  with  increasing  tem¬ 
perature,  especially  for  macropores  (radius  range  from  1  p,m  to 
100  p,m).  They  suggested  that  the  temperature  has  a  small  effect 
on  the  wettability  characteristics  of  MPL  due  to  the  minute  pores 
restricting  the  imbibition  of  water,  therefore  impeding  observation 
of  the  apparent  effects  of  temperature  on  the  pore  characteristics 
of  the  MPL.  However,  the  water  vapor  can  enter  the  minute  pores 
of  the  MPL.  After  the  steam  condensed,  the  liquid  water  remained 
in  the  minute  pores  of  the  MPL.  If  so,  the  temperature  should  also 
have  apparent  effects  on  the  pore  characteristics  of  the  MPL.  Here, 
we  assumed  that  both  the  surface  tension  and  the  contact  angle  of 
the  GDL  as  well  as  the  MPL  decrease  with  increasing  temperature. 
At  353  K,  their  surface  tension  was  0.0625  N  m-1  [2].  The  contact 
angle  was  assumed  to  be  80°  for  the  GDL  and  MPL  surfaces.  Because 
the  values  of  the  surface  tension  and  the  contact  angle  had  already 
incorporated  the  effect  of  temperature  individually,  it  was  not  nec¬ 
essary  to  apply  Eq.  (10)  also.  As  a  special  case,  the  liquid  water 
saturation  was  solved  by  Eq.  (7).  As  shown  in  Fig.  5,  the  satura¬ 
tion  was  significantly  higher  than  that  in  the  GDL  +  MPL  operated 
at  room  temperature.  In  particular,  the  saturation  in  the  MPL  was 
close  to  1,  i.e.,  the  MPL  pores  were  almost  entirely  occupied  by  liq¬ 
uid  water.  In  order  to  maintain  the  capillary  pressure  equilibrium 
at  the  interface  of  the  GDL  and  MPL,  a  high  saturation  is  neces¬ 
sary  in  the  MPL  because  of  its  minute  effective  pore  diameter,  i.e., 
(/C/g)1/2.  Obviously,  serious  flooding  would  occur  in  such  a  case. 
This  is  a  possible  reason  why  Lim  and  Wang  found  that  MPL  has  an 
undesirable  influence  on  the  performance  of  fuel  cells  at  high  cell 
temperatures  and  in  an  over-saturated  environment.  Because  the 
wettability  of  the  GDL  and  MPL  become  hydrophilic,  the  capillary 
pressure  is  negative  at  353  K. 

Finally,  we  would  like  to  emphasize  that  the  phase  change  of 
water  was  neglected  in  this  work  for  the  purpose  of  simplifica¬ 
tion.  We  noted  that  Weber  and  Newman  [28]  have  found  that 
the  evaporation  and  condensation  of  water  caused  a  significant 
movement  of  water,  mainly  through  heat-pipe  effect.  The  phase 
change  of  water  was  also  considered  in  the  models  proposed  by 
Nguyen’s  group  [19,45],  Nam  and  Kaviany  [5],  Mazumder  and 
Cole  [46].  Therefore,  a  coupled  model  including  the  temperature 
influence  on  capillary  pressure  as  well  as  phase  change  is  nec¬ 
essary  for  the  future  research  on  the  water  removal  in  the  fuel 
cells. 


4.  Conclusions 


In  a  one-dimensional  steady-state  model  including  the  effect  of 
temperature  on  the  capillary  pressure,  the  liquid  water  saturation 
and  capillary  pressure  distributions  were  investigated  in  the  cath¬ 
ode  GDL  of  a  PEFC.  The  numerical  results  indicate  that  liquid  water 
saturation  significantly  increases  with  increases  in  the  operating 
temperature  of  the  fuel  cell.  An  elevated  operating  temperature 
has  an  undesirable  influence  on  the  removal  of  liquid  water  inside 
the  GDL.  A  reported  peculiar  phenomenon,  that  the  flooding  of  a 
fuel  cell  operated  under  high  temperature  and  in  an  over-saturated 
environment  is  more  serious  in  a  GDL  combined  with  an  MPL 
than  in  a  GDL  without  an  MPL  is  explained  based  on  the  present 
analysis. 
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